TGFBR1 Variants Can Associate with Non-Syndromic Congenital Heart Disease without Aortopathy.

BACKGROUND: Congenital heart diseases (CHD) are the most common congenital malformations in newborns and remain the leading cause of mortality among infants under one year old. Molecular diagnosis is crucial to evaluate the recurrence risk and to address future prenatal diagnosis. Here, we describe two families with various forms of inherited non-syndromic CHD and the genetic work-up and resultant findings. METHODS: Next-generation sequencing (NGS) was employed in both families to uncover the genetic cause. In addition, we performed functional analysis to investigate the consequences of the identified variants in vitro. RESULTS: NGS identified possible causative variants in both families in the protein kinase domain of the TGFBR1 gene. These variants occurred on the same amino acid, but resulted in differently substituted amino acids (p.R398C/p.R398H). Both variants co-segregate with the disease, are extremely rare or unique, and occur in an evolutionary highly conserved domain of the protein. Furthermore, both variants demonstrated a significantly altered TGFBR1-smad signaling activity. Clinical investigation revealed that none of the carriers had (signs of) aortopathy. CONCLUSION: In conclusion, we describe two families, with various forms of inherited non-syndromic CHD without aortopathies, associated with unique/rare variants in TGFBR1 that display altered TGF-beta signaling. These findings highlight involvement of TGFBR1 in CHD, and warrant consideration of potential causative TGFBR1 variants also in CHD patients without aortopathies.

Identification of novel genetic risk factors of dilated cardiomyopathy: from canine to human.

BACKGROUND: Dilated cardiomyopathy (DCM) is a life-threatening heart disease and a common cause of heart failure due to systolic dysfunction and subsequent left or biventricular dilatation. A significant number of cases have a genetic etiology; however, as a complex disease, the exact genetic risk factors are largely unknown, and many patients remain without a molecular diagnosis. METHODS: We performed GWAS followed by whole-genome, transcriptome, and immunohistochemical analyses in a spontaneously occurring canine model of DCM. Canine gene discovery was followed up in three human DCM cohorts. RESULTS: Our results revealed two independent additive loci associated with the typical DCM phenotype comprising left ventricular systolic dysfunction and dilatation. We highlight two novel candidate genes, RNF207 and PRKAA2, known for their involvement in cardiac action potentials, energy homeostasis, and morphology. We further illustrate the distinct genetic etiologies underlying the typical DCM phenotype and ventricular premature contractions. Finally, we followed up on the canine discoveries in human DCM patients and discovered candidate variants in our two novel genes. CONCLUSIONS: Collectively, our study yields insight into the molecular pathophysiology of DCM and provides a large animal model for preclinical studies.

Genetics of sinoatrial node function and heart rate disorders.

The sinoatrial node (SAN) is the primary pacemaker of the mammalian heart, initiating its electrical activation and ensuring that the heart's functional cardiac output meets physiological demand. SAN dysfunction (SND) can cause complex cardiac arrhythmias that can manifest as severe sinus bradycardia, sinus arrest, chronotropic incompetence and increased susceptibility to atrial fibrillation, among other cardiac conditions. SND has a complex aetiology, with both pre-existing disease and heritable genetic variation predisposing individuals to this pathology. In this Review, we summarize the current understanding of the genetic contributions to SND and the insights that they provide into this disorder's underlying molecular mechanisms. With an improved understanding of these molecular mechanisms, we can improve treatment options for SND patients and develop new therapeutics.

Familial multiple discoid fibromas is linked to a locus on chromosome 5 including the FNIP1 gene.

Previously, we reported a series of families presenting with trichodiscomas, inherited in an autosomal dominant pattern. The phenotype was named familial multiple discoid fibromas (FMDF). The genetic cause of FMDF remained unknown so far. Trichodiscomas are skin lesions previously reported to be part of the same spectrum as the fibrofolliculoma observed in Birt-Hogg-Dubé syndrome (BHD), an inherited disease caused by pathogenic variants in the FLCN gene. Given the clinical and histological differences with BHD and the exclusion of linkage with the FLCN locus, the phenotype was concluded to be distinct from BHD. We performed extensive clinical evaluations and genetic testing in ten families with FMDF. We identified a FNIP1 frameshift variant in nine families and genealogical studies showed common ancestry for eight families. Using whole exome sequencing, we identified six additional rare variants in the haplotype surrounding FNIP1, including a missense variant in the PDGFRB gene that was found to be present in all tested patients with FMDF. Genome-wide linkage analysis showed that the locus on chromosome 5 including FNIP1 was the only region reaching the maximal possible LOD score. We concluded that FMDF is linked to a haplotype on chromosome 5. Additional evaluations in families with FMDF are required to unravel the exact genetic cause underlying the phenotype. When evaluating patients with multiple trichodisomas without a pathogenic variant in the FLCN gene, further genetic testing is warranted and can include analysis of the haplotype on chromosome 5.

Common genetic variants improve risk stratification after the atrial switch operation for transposition of the great arteries.

BACKGROUND: Clinical factors are used to estimate late complication risk in adults after atrial switch operation (AtrSO) for transposition of the great arteries (TGA), but heterogeneity in clinical course remains. We studied whether common genetic variants are associated with outcome and add value to a clinical risk score in TGA-AtrSO patients. METHODS AND RESULTS: This multicenter study followed 133 TGA-AtrSO patients (aged 28 [IQR 24-35] years) for 13 (IQR 9-16) years and examined the association of genome-wide single-nucleotide polymorphisms (SNPs) with a composite endpoint of symptomatic ventricular arrhythmia, heart failure hospitalization, ventricular assist device implantation, heart transplantation, or mortality. Thirty-two patients (24%) reached the endpoint. The genome-wide association study yielded one genome-wide significant (p < 1 × 10-8) locus and 18 suggestive loci (p < 1 × 10-5). A genetic risk score constructed on the basis of independent SNPs with p < 1 × 10-5 was associated with outcome after correction for the clinical risk score (HR = 1.26/point increase [95%CI 1.17-1.35]). Risk stratification improved with a combined risk score (clinical score + genetic score) compared to the clinical score alone (p = 2 × 10-16, C-statistic 0.95 vs 0.85). In 51 patients with a clinical intermediate (5-20%) 5-year risk of events, the combined score reclassified 32 patients to low (<5%) and 5 to high (>20%) risk. Stratified by the combined score, observed 5-year event-free survival was 100%, 79% and 31% for low, intermediate, and high-risk patients, respectively. CONCLUSIONS: Common genetic variants may explain some variation in the clinical course in TGA-AtrSO and improve risk stratification over clinical factors alone, especially in patients at intermediate clinical risk. These findings support the hypothesis that including genetic variants in risk assessment may be beneficial.

PRDM10 directs FLCN expression in a novel disorder overlapping with Birt-Hogg-Dubé syndrome and familial lipomatosis.

Birt-Hogg-Dubé syndrome (BHD) is an autosomal dominant disorder characterized by fibrofolliculomas, pulmonary cysts, pneumothoraces and renal cell carcinomas. Here, we reveal a novel hereditary disorder in a family with skin and mucosal lesions, extensive lipomatosis and renal cell carcinomas. The proband was initially diagnosed with BHD based on the presence of fibrofolliculomas, but no pathogenic germline variant was detected in FLCN, the gene associated with BHD. By whole exome sequencing we identified a heterozygous missense variant (p.(Cys677Tyr)) in a zinc-finger encoding domain of the PRDM10 gene which co-segregated with the phenotype in the family. We show that PRDM10Cys677Tyr loses affinity for a regulatory binding motif in the FLCN promoter, abrogating cellular FLCN mRNA and protein levels. Overexpressing inducible PRDM10Cys677Tyr in renal epithelial cells altered the transcription of multiple genes, showing overlap but also differences with the effects of knocking out FLCN. We propose that PRDM10 controls an extensive gene program and acts as a critical regulator of FLCN gene transcription in human cells. The germline variant PRDM10Cys677Tyr curtails cellular folliculin expression and underlies a distinguishable syndrome characterized by extensive lipomatosis, fibrofolliculomas and renal cell carcinomas.

Biallelic variants in the calpain regulatory subunit CAPNS1 cause pulmonary arterial hypertension.

Purpose: The aim of this study was to identify the monogenic cause of pulmonary arterial hypertension (PAH), a multifactorial and often fatal disease, in 2 unrelated consanguine families. Methods: We performed exome sequencing and validated variant pathogenicity by whole-blood RNA and protein expression analysis in both families. Further RNA sequencing of preserved lung tissue was performed to investigate the consequences on selected genes that are involved in angiogenesis, proliferation, and apoptosis. Results: We identified 2 rare biallelic variants in CAPNS1, encoding the regulatory subunit of calpain. The variants cosegregated with PAH in the families. Both variants lead to loss of function (LoF), which is demonstrated by aberrant splicing resulting in the complete absence of the CAPNS1 protein in affected patients. No other LoF CAPNS1 variant was identified in the genome data of more than 1000 patients with unresolved PAH. Conclusion: The calpain holoenzyme was previously linked to pulmonary vascular development and progression of PAH in patients. We demonstrated that biallelic LoF variants in CAPNS1 can cause idiopathic PAH by the complete absence of CAPNS1 protein. Screening of this gene in patients who are affected by PAH, especially with suspected autosomal recessive inheritance, should be considered.

No prominent role for complement C1-esterase inhibitor in Marfan syndrome mice.

Marfan syndrome (MFS) is a connective tissue disorder causing aortic aneurysm formation. Currently, only prophylactic aortic surgery and blood pressure-lowering drugs are available to reduce the risk of aortic rupture. Upon whole genome sequencing of a Marfan family, we identified a complement gene C1R variant (p.Ser152Leu), which is associated with severe aortic patients. Therefore, we assessed the role of complement activation in MFS aortic tissue. Expression of various complement genes and proteins was detected in human and murine MFS aneurysm tissue, which prompted us to study complement inhibition in MFS mice. Treatment of the Fbn1C1041G/+ MFS mice with human plasma-derived C1-esterase inhibitor Cetor® resulted in reduced complement deposition, decreased macrophage influx in the aorta, and lower circulating TNFα levels. However, in line with previous anti-inflammatory treatments, complement inhibition did not change the aortic dilatation rate in this MFS mouse model. Thus, while complement factors/component 3 activation were detected in human/murine MFS aorta, Cetor® had no effect on aortic dilatation in MFS mice, indicating that complement inhibition is not a suitable treatment strategy in MFS.

Patient-Specific TBX5-G125R Variant Induces Profound Transcriptional Deregulation and Atrial Dysfunction.

BACKGROUND: The pathogenic missense variant p.G125R in TBX5 (T-box transcription factor 5) causes Holt-Oram syndrome (also known as hand-heart syndrome) and early onset of atrial fibrillation. Revealing how an altered key developmental transcription factor modulates cardiac physiology in vivo will provide unique insights into the mechanisms underlying atrial fibrillation in these patients. METHODS: We analyzed ECGs of an extended family pedigree of Holt-Oram syndrome patients. Next, we introduced the TBX5-p.G125R variant in the mouse genome (Tbx5G125R) and performed electrophysiologic analyses (ECG, optical mapping, patch clamp, intracellular calcium measurements), transcriptomics (single-nuclei and tissue RNA sequencing), and epigenetic profiling (assay for transposase-accessible chromatin using sequencing, H3K27ac [histone H3 lysine 27 acetylation] CUT&RUN [cleavage under targets and release under nuclease sequencing]). RESULTS: We discovered high incidence of atrial extra systoles and atrioventricular conduction disturbances in Holt-Oram syndrome patients. Tbx5G125R/+ mice were morphologically unaffected and displayed variable RR intervals, atrial extra systoles, and susceptibility to atrial fibrillation, reminiscent of TBX5-p.G125R patients. Atrial conduction velocity was not affected but systolic and diastolic intracellular calcium concentrations were decreased and action potentials were prolonged in isolated cardiomyocytes of Tbx5G125R/+ mice compared with controls. Transcriptional profiling of atria revealed the most profound transcriptional changes in cardiomyocytes versus other cell types, and identified over a thousand coding and noncoding transcripts that were differentially expressed. Epigenetic profiling uncovered thousands of TBX5-p.G125R-sensitive, putative regulatory elements (including enhancers) that gained accessibility in atrial cardiomyocytes. The majority of sites with increased accessibility were occupied by Tbx5. The small group of sites with reduced accessibility was enriched for DNA-binding motifs of members of the SP (specificity protein) and KLF (Krüppel-like factor) families of transcription factors. These data show that Tbx5-p.G125R induces changes in regulatory element activity, alters transcriptional regulation, and changes cardiomyocyte behavior, possibly caused by altered DNA binding and cooperativity properties. CONCLUSIONS: Our data reveal that a disease-causing missense variant in TBX5 induces profound changes in the atrial transcriptional regulatory network and epigenetic state in vivo, leading to arrhythmia reminiscent of those seen in human TBX5-p.G125R variant carriers.

Correction: Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

[This corrects the article DOI: 10.1371/journal.pgen.1009679.].

Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

Numerous genetic studies have established a role for rare genomic variants in Congenital Heart Disease (CHD) at the copy number variation (CNV) and de novo variant (DNV) level. To identify novel haploinsufficient CHD disease genes, we performed an integrative analysis of CNVs and DNVs identified in probands with CHD including cases with sporadic thoracic aortic aneurysm. We assembled CNV data from 7,958 cases and 14,082 controls and performed a gene-wise analysis of the burden of rare genomic deletions in cases versus controls. In addition, we performed variation rate testing for DNVs identified in 2,489 parent-offspring trios. Our analysis revealed 21 genes which were significantly affected by rare CNVs and/or DNVs in probands. Fourteen of these genes have previously been associated with CHD while the remaining genes (FEZ1, MYO16, ARID1B, NALCN, WAC, KDM5B and WHSC1) have only been associated in small cases series or show new associations with CHD. In addition, a systems level analysis revealed affected protein-protein interaction networks involved in Notch signaling pathway, heart morphogenesis, DNA repair and cilia/centrosome function. Taken together, this approach highlights the importance of re-analyzing existing datasets to strengthen disease association and identify novel disease genes and pathways.

Terminal osseous dysplasia with pigmentary defects and cardiomyopathy caused by a novel FLNA variant.

Terminal osseous dysplasia with pigmentary defects (TODPD), also known as digitocutaneous dysplasia, is one of the X-linked filaminopathies caused by a variety of FLNA-variants. TODPD is characterized by skeletal defects, skin fibromata and dysmorphic facial features. So far, only a single recurrent variant (c.5217G>A;p.Val1724_Thr1739del) in FLNA has found to be responsible for TODPD. We identified a novel c.5217+5G>C variant in FLNA in a female proband with skeletal defects, skin fibromata, interstitial lung disease, epilepsy, and restrictive cardiomyopathy. This variant causes mis-splicing of exon 31 predicting the production of a FLNA-protein with an in-frame-deletion of 16 residues identical to the miss-splicing-effect of the recurrent TODPD c.5217G>A variant. This mis-spliced transcript was explicitly detected in heart tissue, but was absent from blood, skin, and lung. X-inactivation analyses showed extreme skewing with almost complete inactivation of the mutated allele (>90%) in these tissues, except for heart. The mother of the proband, who also has fibromata and skeletal abnormalities, is also carrier of the FLNA-variant and was diagnosed with noncompaction cardiomyopathy after cardiac screening. No other relevant variants in cardiomyopathy-related genes were found. Here we describe a novel variant in FLNA (c.5217+5G>C) as the second pathogenic variant responsible for TODPD. Cardiomyopathy has not been described as a phenotypic feature of TODPD before.

Rare variants in KDR, encoding VEGF Receptor 2, are associated with tetralogy of Fallot.

PURPOSE: Rare genetic variants in KDR, encoding the vascular endothelial growth factor receptor 2 (VEGFR2), have been reported in patients with tetralogy of Fallot (TOF). However, their role in disease causality and pathogenesis remains unclear. METHODS: We conducted exome sequencing in a familial case of TOF and large-scale genetic studies, including burden testing, in >1,500 patients with TOF. We studied gene-targeted mice and conducted cell-based assays to explore the role of KDR genetic variation in the etiology of TOF. RESULTS: Exome sequencing in a family with two siblings affected by TOF revealed biallelic missense variants in KDR. Studies in knock-in mice and in HEK 293T cells identified embryonic lethality for one variant when occurring in the homozygous state, and a significantly reduced VEGFR2 phosphorylation for both variants. Rare variant burden analysis conducted in a set of 1,569 patients of European descent with TOF identified a 46-fold enrichment of protein-truncating variants (PTVs) in TOF cases compared to controls (P = 7 × 10-11). CONCLUSION: Rare KDR variants, in particular PTVs, strongly associate with TOF, likely in the setting of different inheritance patterns. Supported by genetic and in vivo and in vitro functional analysis, we propose loss-of-function of VEGFR2 as one of the mechanisms involved in the pathogenesis of TOF.

Biallelic loss-of-function variants in PLD1 cause congenital right-sided cardiac valve defects and neonatal cardiomyopathy.

Congenital heart disease is the most common type of birth defect, accounting for one-third of all congenital anomalies. Using whole-exome sequencing of 2718 patients with congenital heart disease and a search in GeneMatcher, we identified 30 patients from 21 unrelated families of different ancestries with biallelic phospholipase D1 (PLD1) variants who presented predominantly with congenital cardiac valve defects. We also associated recessive PLD1 variants with isolated neonatal cardiomyopathy. Furthermore, we established that p.I668F is a founder variant among Ashkenazi Jews (allele frequency of ~2%) and describe the phenotypic spectrum of PLD1-associated congenital heart defects. PLD1 missense variants were overrepresented in regions of the protein critical for catalytic activity, and, correspondingly, we observed a strong reduction in enzymatic activity for most of the mutant proteins in an enzymatic assay. Finally, we demonstrate that PLD1 inhibition decreased endothelial-mesenchymal transition, an established pivotal early step in valvulogenesis. In conclusion, our study provides a more detailed understanding of disease mechanisms and phenotypic expression associated with PLD1 loss of function.

Germline variants in HEY2 functional domains lead to congenital heart defects and thoracic aortic aneurysms.

PURPOSE: In this study we aimed to establish the genetic cause of a myriad of cardiovascular defects prevalent in individuals from a genetically isolated population, who were found to share a common ancestor in 1728. METHODS: Trio genome sequencing was carried out in an index patient with critical congenital heart disease (CHD); family members had either exome or Sanger sequencing. To confirm enrichment, we performed a gene-based association test and meta-analysis in two independent validation cohorts: one with 2685 CHD cases versus 4370 . These controls were also ancestry-matched (same as FTAA controls), and the other with 326 cases with familial thoracic aortic aneurysms (FTAA) and dissections versus 570 ancestry-matched controls. Functional consequences of identified variants were evaluated using expression studies. RESULTS: We identified a loss-of-function variant in the Notch target transcription factor-encoding gene HEY2. The homozygous state (n = 3) causes life-threatening congenital heart defects, while 80% of heterozygous carriers (n = 20) had cardiovascular defects, mainly CHD and FTAA of the ascending aorta. We confirm enrichment of rare risk variants in HEY2 functional domains after meta-analysis (MetaSKAT p = 0.018). Furthermore, we show that several identified variants lead to dysregulation of repression by HEY2. CONCLUSION: A homozygous germline loss-of-function variant in HEY2 leads to critical CHD. The majority of heterozygotes show a myriad of cardiovascular defects.

Expanding the phenotype of biallelic RNPC3 variants associated with growth hormone deficiency.

Pathogenic variants in components of the minor spliceosome have been associated with several human diseases. Recently, it was reported that biallelic RNPC3 variants lead to severe isolated growth hormone deficiency and pituitary hypoplasia. The RNPC3 gene codes for the U11/U12-65K protein, a component of the minor spliceosome. The minor spliceosome plays a role in the splicing of minor (U12-type) introns, which are present in ~700-800 genes in humans and represent about 0.35% of all introns. Here, we report a second family with biallelic RNPC3 variants in three siblings with a growth hormone deficiency, central congenital hypothyroidism, congenital cataract, developmental delay/intellectual deficiency and delayed puberty. These cases further confirm the association between biallelic RNPC3 variants and severe postnatal growth retardation due to growth hormone deficiency. Furthermore, these cases show that the phenotype of this minor spliceosome-related disease might be broader than previously described.

Loss-of-function variants in myocardin cause congenital megabladder in humans and mice.

Myocardin (MYOCD) is the founding member of a class of transcriptional coactivators that bind the serum-response factor to activate gene expression programs critical in smooth muscle (SM) and cardiac muscle development. Insights into the molecular functions of MYOCD have been obtained from cell culture studies, and to date, knowledge about in vivo roles of MYOCD comes exclusively from experimental animals. Here, we defined an often lethal congenital human disease associated with inheritance of pathogenic MYOCD variants. This disease manifested as a massively dilated urinary bladder, or megabladder, with disrupted SM in its wall. We provided evidence that monoallelic loss-of-function variants in MYOCD caused congenital megabladder in males only, whereas biallelic variants were associated with disease in both sexes, with a phenotype additionally involving the cardiovascular system. These results were supported by cosegregation of MYOCD variants with the phenotype in 4 unrelated families by in vitro transactivation studies in which pathogenic variants resulted in abrogated SM gene expression and by the finding of megabladder in 2 distinct mouse models with reduced Myocd activity. In conclusion, we have demonstrated that variants in MYOCD result in human disease, and the collective findings highlight a vital role for MYOCD in mammalian organogenesis.

GATA6 mutations: Characterization of two novel patients and a comprehensive overview of the GATA6 genotypic and phenotypic spectrum.

The first human mutations in GATA6 were described in a cohort of patients with persistent truncus arteriosus, and the phenotypic spectrum has expanded since then. This study underscores the broad phenotypic spectrum by presenting two patients with de novo GATA6 mutations, both exhibiting complex cardiac defects, pancreatic, and other abnormalities. Furthermore, we provided a detailed overview of all published human genetic variation in/near GATA6 published to date and the associated phenotypes (n = 78). We conclude that the most common phenotypes associated with a mutation in GATA6 were structural cardiac and pancreatic abnormalities, with a penetrance of 87 and 60%, respectively. Other common malformations were gallbladder agenesis, congenital diaphragmatic hernia, and neurocognitive abnormalities, mostly developmental delay. Fifty-eight percent of the mutations were de novo, and these patients more often had an anomaly of intracardiac connections, an anomaly of the great arteries, and hypothyroidism, compared with those with inherited mutations. Functional studies mostly support loss-of-function as the pathophysiological mechanism. In conclusion, GATA6 mutations give a wide range of phenotypic defects, most frequently malformations of the heart and pancreas. This highlights the importance of detailed clinical evaluation of identified carriers to evaluate their full phenotypic spectrum.